Brown-HET-1586 



January 18, 2011 

Dynamics of the chiral phase transition from AdS/CFT duality 

Gerald Guralnik/'Q Zachary Guralnik/'Q and Cengiz Pehlevan 2, 

1 Department of Physics, Brown University, Providence, RI 02906 USA 
2 Edmond and Lily Center for Brain Sciences, 
The Interdisciplinary Center For Neural Computation, 
The Hebrew University of Jerusalem, Jerusalem, 91904 Israel 

Abstract 

We use Lorentzian signature AdS / CFT duality to study a first order phase transition in strongly 
coupled gauge theories which is akin to the chiral phase transition in QCD. We discuss the relation 
between the latent heat and the energy (suitably defined) of the component of a D-brane which 
lies behind the horizon at the critical temperature. A numerical simulation of a dynamical phase 
transition in an expanding, cooling Quark-Gluon plasma produced in a relativistic collision is 
carried out. 
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I. INTRODUCTION 



At present, dual string theory descriptions of a number of large N gauge theories are 



known and facilitate computations at large 't Hooft coupling via AdS/CFT duality 
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431 ] . Some of these are distant 1 cousins of QCD, and capture some of the qualitative 



features of strong coupling phenomena in QCD. Particular attention has been given to the 
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phase structure of these theories i 

in heavy ion collisions si, [2^]. A prototype of the chiral phase transition in QCD is a 



181 ] and the behavior under conditions arising 



first order phase transition in a strong 



theory, discovered in 



flflflfl 
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coupled large N N — 2 supersymmetric gauge 
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using the dual AdS description. In the 



Euclidean description, this transition is realized as a change in topology of space filling D7- 
branes in an asymptotically AdS-black hole background, and corresponds to a jump in a 
ipif) condensate. A background dual to a theory more like QCD is described in 41], in which 
case the analogous change in topology corresponds to a chiral symmetry breaking transition. 

In this article, we revisit the phase transition in the AdS dual of the prototype M = 2 
supersymmetric model from a Lorentzian signature perspective, with an eye toward study- 
ing time dependent non-equilibrium processes. Non-equilibrium processes, such as the ex- 
pansion and cooling of a quark-gluon fireball produced in a heavy ion collision, are not 
amenable to lattice gauge theory methods, but may be approached via a Lorentzian signa- 
ture AdS/CFT duality. Furthermore, there are very interesting questions about the inter- 
pretation of thermodynamic quantities in the dual gravitational description which are better 
asked in Lorentzian signature. While quantities such as entropy and latent heat can be com- 
puted semi-classically from the Euclidean gravity/D-brane action, they lack a satisfactory 
physical interpretation 2 . 

In Euclidean signature, the boundary of the AdS-Schwarzchild black hole has topology 
S* 5 <S> S 1 (g) R 3 which, in the M = 2 model, is wrapped by D7-branes having the boundary 
topology S 3 £§> S 1 <E> -R 3 - The S 1 of the space-time contracts to zero size in the interior, 
at a point corresponding to the horizon after continuation to Lorentzian signature. In the 
low temperature phase the D7-brane ends smoothly in the interior by contraction of the S 3 



1 Or close, depending on which question is being asked 

2 The Euclidean results correspond to an assumed saddle point approximation to integrals over degrees of 
freedom (e.g. string fields) which are not specified or even completely understood 
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wrapping the space-time S 5 . In the high temperature phase the D7-brane ends smoothly 
at the same point where the S l of the Euclidean space-time, which the D7-brane wraps, 
contracts to zero size. 

Upon continuation to Lorentzian signature, the D7-brane in the low temperature phase 
ends before reaching the horizon, while the D7-brane extends through the horizon in the high 
temperature phase. We will solve the static equations of motion for the D7-brane in the in- 
falling case in Section UTTl including the region behind the horizon, where it can be shown that 
the D7-brane ends at a conical singularity. Since the AdS radius becomes time-like inside 
the horizon, the conical singularity can be viewed as one-half of an annihilation diagram: 
the D7-brane ends before reaching the black-hole singularity by annihilation into closed- 
string modes. We show that the conical singularity corresponds to a spatial S 1 collapsing 
at the speed of light. In passing dynamically from the high temperature phase to the low 
temperature phase, the D7-brane can be said to "pull out" of the black hole, although what 
in fact is happening is that the component of the D7-brane interior to the horizon vanishes 
by annihilation into gravitons, etc. 

While for many questions the component of the D7-brane interior to the horizon is irrel- 
evant, it is presumably related to the existence of latent heat in the AdS-description of the 
phase transition. In the low temperature phase the free energy, obtained from the Euclidean 
action, can be shown to be equivalent to a suitably defined energy of the D7-brane in the 
Lorentzian signature description, while in the high temperature phase the free energy corre- 
sponds to the energy of the component of the D7-brane lying outside the horizon. The free 
energy is a bound on the energy available for work, which in this case corresponds to flux 
of energy across the AdS boundary. This flux is bounded by the energy of the component 
of the D7-brane lying outside the horizon. The flux of energy across the AdS boundary is 
defined only after a suitable renormalization, discussed in Section IIV1 In a static AdS-black 
hole background, and neglecting back-reaction, there is conserved energy-momentum tensor 
associated with time translation invariance 3 , such that passing between two configurations 
with the same free energy requires flux across the AdS boundary match flux across the 
horizon. The energy of the component of the D7-brane inside the horizon in the high tem- 
perature phase would therefore seem to be related to the amount of energy required for the 

3 Inside the horizon, the variable in which the metric is translation invariant becomes space-like. 
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transition from low temperature (non-inf ailing) to the high temperature (infalling) phase, 
i.e. the latent heat. We evaluate the energy of the component of the D7-brane interior to 
the horizon at the critical temperature in section HV C\ finding a result which is about 1/8 
of the latent heat determined by the discontinuity in the Tc^J 7 across the phase transition. 
We do not yet have a satisfactory understanding of the discrepancy between the interior 
energy and the latent heat, although it is presumably related to the neglected back reaction. 
With an eye towards simulating a strongly coupled expanding-cooling plasma, we con- 



struct a metric dual to an expanding cooling p 
constructed previously; the initial attempts 



asma in section |VJ Similar metrics have been 



25| having a singularity where there should 



have been a horizon. This problem was rectified in [23|, |29|, using Eddington-Finkelstein 
type coordinates. We find it much more convenient to use a Painleve-Gullstrand or 'river 
model' type metric. In section IVIt we describe our initial efforts to simulate the dynamics 
of a D7-brane passing through the first order phase transition in this background. 

The work presented here appeared before in one of the authors' PhD thesis 40] . A related 
paper [yj] appeared as we were finalizing this work. 



II. EQUILIBRIUM THERMODYNAMICS IN EUCLIDEAN SIGNATURE 



We consider the M = 2 super Yang-Mills theory obtained by coupling an M = 4 multiplet 
to Nf hypermultiplets in the fundamental representation of the gauge group SU(N C ). This 
theory is conformal in the limit N c — > oo with Nf fixed 4 . For large t'Hooft coupling A = 
g 2 N c » 1, the theory is dual to weakly coupled supergravity in an Anti-deSitter background 
with space-filling D7-branes 26|. The phase diagram and spectroscopy of this theory has 
been studied in great detail using AdS / CFT duality. Of particular note is a first order phase 



transition as the ratio of temperature to c 



uark mass is varied. This transition is akin to the 



chiral phase transition in QCD [3, l6|, l8|, 130(] . 

Equilibrium thermodynamic quantities in this theory can be computed from the dual 
string theory description of this theory, which involves a background with a Euclidean black 



For finite N c , this theory is not asymptotically free and requires an ultraviolet completion. 
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hole in AdS 5 x S 5 space time, 

ds 2 = 1 ((1 - b A z A )dt 2 + + dx 2 ) + dQ 2 , (1) 

where it is convenient to write the five-sphere metric as 

dVtl = d6 2 + cos 2 6{d<f) 2 + cos 2 <j)dQ 2 3 ) . (2) 

Throughout this paper, we will set the AdS curvature, R, to 1. This space-time is defined 
for z < 1/b and is smooth and complete if the Euclidean time t is compactified on a circle 
of radius ir/b, corresponding to the inverse temperature. 

The inclusion of the Nf hypermultiplets (quarks) corresponds to the addition of Nf D7- 
branes to this background, embedded on a surface 

= 0, = 6(z). (3) 

The induced metric on the D7-brane is 

ds 2 = -i ( (1 - b 4 z 4 )dt 2 + \r-rdz 2 + dx 2 ] + O'izfdz 2 + cos 2 OdSj , (4) 

z 2 \ 1 — 4 Z 4 / 



giving the D7-brane action 



I m = -Mp I dz^^-^l + z 2 (l-Wz*)9' 2 . (5) 



horizon Z 




with Af = NfTwifls (J d 3 x), where Td7 is the brane tension and Q3 is the volume of the 
three-sphere. The action ((Sj) is divergent as the boundary of integration z — > and needs 
to be renormalized, which will be discussed in more detail below. This gives rise to the 
equations of motion 

=^ 2 (-l + b A z 4 )6" + z(3 + b A z 4 )6' + 2z 3 (l - 6V)(2 - feV)# /3 

+ 3tan# (-1 + z 2 (-l + b A Z y 2 ) . (6) 

Near the AdS boundary, z — > 0, the solutions have the asymptotic behavior 

6{z) ~ B x bz + 9 3 b 3 z 3 + • ■ ■ , (7) 

where the parameters B\ 3 determine the quark mass M and chiral condensate C. Specifically 



m = l -V\Te u (8) 

C = (qq) = -lN f N c T 3 VX f-29 3 + |) . (9) 



The topology of a constant z slice of the D7-brane at the AdS boundary z — > is t&^S 1 ® 
S 3 . The D7-brane can end at z end < 1/6 if the S 3 contracts to zero size, with #(z en d) = 7r/2, 
and ^'(^end) = oo so that the end is smooth rather than a conical singularity. If on the other 
hand z en & = 1/6, the S 1 of the space-time in which the D7-brane is embedded contracts 
to zero and the boundary condition 9'(z end ) = |tan^(z end ) follows from the equations of 
motion ([6]). Solving the D7-brane equations of motion at fixed temperature, subject to the 
boundary conditions at z end , yields a curve in the C-M plane parameterized by z end for 
D7-branes which do not extend to the horizon and #(zhorizon) for D7-branes which end at 
the horizon (see figured]). For these embeddings, the free energy F is given by film, and 
the free energy density T = Fj J d 3 x satisfies 

dT = CdM (10) 

at constant temperature. In a certain range of values for M there are solutions with several 
values of C, indicating a first order phase transition. The transition occurs at a critical 
value of the only dimensionless parameter T/M. At a given temperature, the critical mass 
can be determined by Maxwell's equal area construction |2j. It lies the a point M = M c 
where there are three solutions, such the integral over the curve C(z cnd ), M(z en< i) connecting 
them gives J CdM = 0. The physical solutions at the critical point are the outer two which 
for which the free energy is convex, -jrh = ^ > 0, and takes the same value via the 
equal area rule. Physical transitions from one phase to the other do not follow the curve 
C(Zead), ^(^end), but occur by bubble nucleation or supercooling/heating. The solutions we 
describe subsequently focus on supercooling/heating which is more numerically tractable. 

III. STATIC D7-BRANE EMBEDDING IN LORENTZIAN SIGNATURE 

Before considering D7-branes embedded in the time dependent background (1351) . we will 
discuss some of the properties of the static equilibrium solution in Lorentzian signature. 
While the static embedding outside the horizon is the same in either Euclidean or Lorentzian 
signature, it is interesting to consider the behavior of the D7-branes behind the horizon, a 
region which exists only in Lorentzian signature. For reasons to be discussed later, the 
behavior of the D7-brane behind the horizon is related to the latent heat. 
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FIG. 1: On top, D7-brane embeddings in the Euclidean black hole background is shown, b is set 
to 1. Blue (solid) lines plot embeddings that extend to the horizon and red (dashed) lines plot 
embeddings that do not extend to the horizon. Below, quark condensate as a function of M/T 
that arises from the D7-brane embeddings is shown. The figure on bottom right zooms into the 
multivalued region of the figure on bottom left. The green (vertical) line shows the critical mass 
at which the condensates vev jumps discontinously, at which t^t^ = 0.9234. This value was found 
by equating the area between the green (vertical) line and the CM curve on both sides of the green 
(vertical) line [3]. 

The Lorentzian signature AdS-Schwarzchild black hole metric is 

ds 2 = i ^-(1 - b A z A )dt 2 + 1 _\ 4z4 dz 2 + dx 2 ^j + dn 2 5 . (11) 
A redefinition of the time coordinate, 



gives a metric of the Painleve-Gullstrand, or 'river-model' form 

ds 2 = \ {-dtl + (dz - b 2 z 2 dt r ) 2 + dx 2 ) + dd 2 + cos 2 6(d<p 2 + cos 2 <pdSl) , (13) 

which has the advantage that there is no coordinate singularity at the horizon. Moreover, 
time dependence can later be introduced in the parameter b without generating a curvature 
singularity. Since only the time coordinate has been redefined, the equations for 6(z) for a 
static D7-brane embedding are identical to those in equation (J6j). The only difference from 
the solutions in the Euclidean background ([T]) is that the solutions which reach the horizon 
may be continued behind the horizon to z > 1/6. This is done by enforcing continuity across 
the horizon, see Figure [2j Numerically, one finds that the infalling solutions reach 9 = n/2 
at a finite value of z > 1. Thus, the infalling D7-branes end like the non-infalling D7-branes, 
via an S 3 contracting to zero size. However the ending is not smooth due to causality. Inside 
the horizon z is timelike, so the boundary condition required for smoothness, 9'(z en d) = °°> 
would correspond to an S 3 collapsing faster than the speed of light. Upon inspection of (El), 
the correct boundary condition for infalling solutions is seen to be 

-l + 4^-1 + b'zl A )e\z end ) 2 = o , (14) 

which corresponds to an S 3 collapsing at the speed of light. A suitably defined energy 
associated with the cone lying behind the horizon would seem to be a natural guess for the 
latent heat 5 . We will explore this proposal more deeply in the subsequent discussion. 

IV. REMARKS ON FREE ENERGY AND THE STRESS-ENERGY TENSOR 

In this section we calculate the renormalized stress-energy tensor of the D7-brane. We 
will see that the free energy is equivalent to energy of the component of the static D7-brane 
lying outside the horizon. The free energy represents a bound on mechanical work which can 
be obtained from a thermodynamic system at constant temperature. In the present case, the 
energy external to the horizon represents a bound on energy flux at the AdS boundary. We 
will also compute the component of the D7-brane stress tensor determining the energy flux 
at the AdS-boundary boundary and horizon-crossing. For quasi-static time evolution, we 

5 'Latent' means hidden, which in this case has a very literal realization as hidden behind the horizon 
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FIG. 2: Static D7-brane embedding in Lorentzian signature, b is set to 1. 



will see that the flux at the AdS boundary is CM, consistent with ffHJj) . Although J CdM 
vanishes along the curve of (not always stable) solution connecting the two phases which 
coexist at the critical temperature, the physical passage between phases does not occur this 
way, nor is it quasi-static, such that the net flux at infinity may be non-zero. In a transition 
between these phases which satisfies the equations of motion, there is a lower bound on the 
D7 energy flux at infinity determined by the energy of the cone interior to the horizon in 
the infalling phase. A lower bound on flux at the AdS boundary seems clearly related to 
the existence of latent heat. However we will see that the cone energy accounts for only a 
fraction of the latent heat. 
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A. Renormalized action 



The Euclidean action of the D7-brane, given by (|3j) is divergent as z — > 0, but rendered 



finite by holographic renormalization 



IT. 



42] . Necessary counterterms to renormalize 



the action were given in [27] . One introduces a finite cut-off at z = e slice and adds the 
counterterms made of curvature invariants at the z = e slice. The authors of 27] calculated 



these counterterms in Fefferman-Graham coordinates. The counterterms are given by 



U 

u 

L F 



□ 7 + -It 



6{x,e) 



(15) 



jij is the induced brane metric on the cut-off slice (at coordinate e in the radial direction) 



in these coordinates. We denote the renormalized action by J ren . In 19] . the D7 brane 
action renormalization was done for a time- dependent embedding. The free energy is given 
by F = Iren/P- 111 Lorentzian signature, this turns out to be equivalent to energy (suitably 
defined) of the component of the D7-brane outside the horizon. 



B. Stress energy tensor 

The Hilbert stress energy tensor is 



T AB (x) = 2 * Im . (16) 
V^gog A B 

Given a killing vector £ , one can define conserved currents 

d A (£ B V=9~T A B ) = 0. (17) 
We consider a time dependent embedding, 6 = 9(z,t), on the Lorentzian AdS-black hole 



background, 



ds 2 = 4 f"(l - b A z A )dt 2 + - \ A A dz 2 + dx 2 ) + dQ 2 5 , (18) 
z 2 \ 1 — o 4 ^ 4 
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and the Killing vector £ a 8a = Jp We will call these coordinates " Schwarzschild coordinates" 
in the following. Then, in the absence of counterterms, 

\Z=9 T \ r/,Nr/„ ~/ , r, COS 3 9 1 + Z 2 (l - 



NfT D7 ' Z 5 A , „ 2 ft l.4~4\n, 2 ^ 



*f ±D7 z Jl + z*{l-b*z*)& 



l-b 4 z 4 



- v ^ * = 5 V 5(0 - *)) det fi 3 — , 1 J . • (19) 

N f Tm Z ^l + z\l-Wz±)9>*--^- A 

Here det Q3 refers to the determinant of the three-sphere part of the metric (jSJ) and B(z,t) 

is the classical trajectory. As expected, for static configurations \J—g T z t vanishes. One can 

then define a conserved energy associated with the probe D7-brane, 

/",,,, , . Kr f , cos 3 9 l + z 2 (l-6V)G' 2 , . 

U D7 = - / dzd 3 a;d . . . ^/~ z gT\ = M \ dz = v ; =. (20) 



z " -/l + z 2 (l-6 4 ^)0' 2 -^ 



Thus, when the embedding is static, the conserved energy of the part of the D7 brane outside 
the horizon is formally equal to the brane's free energy (see (JSJ)). 

The stress-energy tensor as defined above is divergent as z — > and leads to divergent 
energies. Including the counterterms (fT5|) in the renormalized action yields the canonical 
stress energy tensor components 6 

V ,., ,™ n , xx , „ cos 3 9 l+z 2 (l-b 4 z 4 )9' 2 



5(ip)8(9 -B(z,t)) det fi 3 



N f T ° 7 ^ Jl + z 2 (l-b*z 



y5 



A\a>2 



z 2 e 2 

l-b 4 z 4 



lim <^ 5(z - e)5(m(9 - Bit, z)W(. . .)£l 3 y/=7 



4 2 12 2 ' v ; 



(21) 

where 7^ is the induced metric on 2; = e slice. Integrating (12T]) over 2; yields a finite 
Hamiltonian, H = f dzT*. There is a conservation law 

d t r z - d z r t z , (22) 

which can be used to renormalize T t z as z — > 

7? , | r , AP //, ~, xx, ^ cos 3 # -z 2 (l-6V)^ 



lim - * = lim < 5(ip)5(9 - B(z, t)) det fi. 



3- 



5 



-5(V)<K0 - e))5 3 (. • .)#UV=7 



3 



(23) 



6 The Hilbert stress tensor computed with the renormalized action does not yield finite energy unless the 
configuration is static. 
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and to show that the flow of energy at the AdS boundary is 

h 2 ... 1 

(24) 



- ton = 6(^)5(6 - Q(t, e))5 3 (. . .)0 3 



ft 4 • • b 2 ■ ■■' 

— 010? - 26 4 0!0 3 010! 

3 1 2 



In a quasi-static process, keeping only terms which are first order in time derivatives, 
Q23D and (ED implies 

H = J CM, (25) 

Note that ( I24p assumes that there can only be flow from the AdS boundary and not the 
endpoint of the D7-brane. This is true if the D7-brane ends (smoothly) before reaching 
the horizon. However there may be flow across the singular endpoint of a D7-brane ending 
inside the horizon. Note that the energy flux across the horizon is 0(6 2 ), and so vanishes in 
quasi-static evolution. This is consistent with dF = CdM where F, the free energy, is the 
component of the energy outside the horizon f° 77*. 

^horizon 

C. Latent Heat 

The two phases which coexist at the critical temperature have the same free energy, 
although the D7-brane Hamiltonian differs by the energy of the cone behind the horizon in 
the infalling phase. Passing dynamically from one phase to the other is not a quasistatic 
process; the energy of the cone 

P ^horizon 

£ C onc= / dzT t \ (26) 

is a bound on the flux of D7-brane energy across the horizon in passing from the non-infalling 
phase to the infalling phase. Given the equivalence of the free energy on either side of the 
transition, this is also a bound on the flux across the AdS boundary. As a bound on the 
energy which must be added to the system in a phase change at the critical temperature, and 
on TAS corresponding to the change in black hole mass, this would seem to be equivalent 
to the latent heat. However this is not the most stringent bound. The latent heat can 
be reliably computed from the discontinuity in the derivative of the free free energy with 
respect to temperature, 

Q latent = TAS = -TA (|0 (27) 
= MAC = 0.0011 XT 4 NfN c . (28) 
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This result is about eight times larger than (126]) . which is computed numerically for the 
infalling solution at the critical temperature. The discrepancy between fl26|) and (1271) is 
presumably related to the fact that back- reaction has not been taken into account. Due to 
the equations of motion, the effect of the back-reaction on the free energy is of smaller order 
in the 1/N C expansion than the effect of the back reaction on the internal energy U = F+TS. 
The effect of back-reaction on the latter is of the same order, NfN c , as the D7-brane energy 



281 ] . The latent heat may be reliably calculated from ( 127)) . while ( 126)) captures only a partial 



contribution. 



V. DYNAMICS OF AN EXPANDING PLASMA 



One of advantages of AdS/CFT duality over lattice gauge theory is the possibility of 
simulating time dependent, non-equilibrium processes in a strongly gauge theory. In the 
subsequent discussion we describe initial efforts at a numerical simulation of the dynamical 
passage through the phase transition in an expanding cooling plasma. Much of this work 
has already appeared in Q]. As this article was being completed, a closely related article 



appeared 



An approximate Lorentzian signature solution of Einstein's equations dual to a boost 
invariant expanding cooling Af = A Yang-Mills plasma, akin to that arising in heavy ion 
collisions, was given in [25|. In Fefferman Graham coordinates, the metric is 

2 



ds 2 = ^ 



1 ep 

1 3 f 4 / 3 



1 _l_ e z 

1 ~r 3 ?*73 



df 2 + 1 



e 2_ 
3 f 4 / 3 



r drj + dx A 



dz 2 

+ 4r , (29) 



where r, rj, x± are the coordinates of the the Yang-Mills theory and of the AdS boundary 
at z — > 0. These are the natural coordinates for a relativistic heavy ion collision where, for 



\\d.jJ^s (the rapidity), and x± = x Y ,x 2 . 



a collision along the x axis, t — i >. ■< ! • if - ... 
The asymptotic z — >■ expansion of the metric yields the expectation value of the Yang- 
Mills theory stress energy tensor, via g^ v = gffl + z 4 < T^ v > + • • • . The stress energy 
tensor derived from (129)) corresponds to a relativistic perfect fluid with boost invariant 



initial conditions 



251 ]. with energy density 



e r 



e 

4/3 ' 



T 



(30) 
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The metric (129]) above is the leading term in a late time expansion in terms of a scaling 
variable s = zf' 1 ^ 3 , such that = g$(s) + 0(f~ 2 ^ 3 ). The dynamics of non-infalling 
D7-branes in this background was studied in 19] . at temperatures well above the phase 
transition. However ( 129"j) is not well suited for studying the non-equilibrium dynamics of the 
phase transition, since it does not give a valid description of the region of space-time in the 
neighborhood of the horizon. This problem was resolved in {23, 29], in which a solution was 
given in Eddington-Finkelstein coordinates. However, we have found yet another form of 
the solution, involving a metric of the Painleve-Gullstrand from, to be far more convenient. 



Proposed Background 



In this section we introduce our proposed background using phenomenological arguments. 
Later on we compare to the metric of Janik and Peschanski, (129]) . and show that both metrics 
are related by a coordinate transformation. Before making our educated guess, let us note 
that the solution (I2"9l closely resembles an AdS-black hole in Fefferman- Graham coordinates, 



z 2 



" b Z ,A) -dt 2 + {l + b i z i /A)dx 2 



+ 



dz 



~2 



-9 ' 
Z l 



(31) 



1 + b i z i /4 

where the replacements t — > f and dx 2 — > f 2 dr] 2 + dx]_ are made to realize boost invariance 
in an expanding plasma, and h is a function of r, 



b(f) 



4eo 
3 



1/4 



-1/3 



(32) 



We require a Lorentzian signature solution which is valid from the boundary through the 
horizon. Therefore, as a starting point consider the following metric of a static AdS black 
hole, 



ds 2 = — (-dt 2 + (dz - b 2 z 2 dt r ) 2 + dx 2 ) 



(33) 



where (133]) is related to the AdS-Schwarzchild metric (TTTj) by the coordinate transformation 

b 2 z' 2 



t - 



rdz' . 



(34) 



1 - b A z' A 

These coordinates are similar to the "river model" or Gullstrand-Painleve coordinates of a 
Schwarzchild black hole in an asymptotically Minkowski space 



21 



To obtain a result 
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consistent with an expanding plasma akin to one produced in a heavy ion collusion, we 
replace the metric flMl) with a boost invariant version and allow b to depend on r; 



ds 2 = \ {-dr 2 + (dz - b(r) 2 z 2 dr) 2 + r 2 dif + dx\) . (35) 

The geometry (|35|) looks like a AdS black hole solution with horizon moving in the bulk 
according to Zh = l/6(r). With an abuse of terminology, we call this surface a horizon, even 
though in the r dependent case it is non-trivial to determine the existence and location of 
the horizon. For slowly varying b(r), the volume form on the horizon is 

A = b(r) 3 rdr] A dx 2 A dx 3 . (36) 

The volume corresponds to an entropy per unit rapidity per unit transverse area S(t) = 
6(r) 3 r, which must increase or remain constant with increasing r. If the expansion is 
isentropic, then b(r) ~ r -1 / 3 . 

With b given by ([3"2l . it is not hard to see that the following change of coordinates 

r = f(l + ^), (37) 



s/TTWTi 

where s = (^) 1//4 zf _1 / 3 and 



T 



4e V 1/4 f s s' 2 



T2 (s) = - — ^ ds' ^ (38) 

K} V 3 J J (l + s'Y4) 1/2 (l-s'Y4) 1 ; 

transforms between our proposal and ( 1291) . to leading order in the f~ 2//3 expansion of 25]. 

The advantage of the metric fl3"5"|) which we propose is that it extends through the horizon 
can therefore be used to study chiral dynamics near the phase transition. Other proposed 



metrics written in Eddington-Finkelstein coordinates [23j, |29| also extend across the horizon, 
however our metric has the convenient property that it is manifestly AdS as z — > 0, which 
is makes it easier to apply the AdS/CFT prescription. 

While we did not obtain ( 1351) by solving Einstein's equations in a controlled expansion, 
although that is presumably possible 7 , we adopt it as the simplest metric with which to study 
non-equilibrium properties of the chiral phase transition in an expanding cooling plasma. 
In particular, we will use this background to numerically simulate the evolution of a probe 
D7-brane 'pulling out' of the black hole. 



7 A possible starting point would be to seek solutions of Einstein's equations of motion with a metric ansatz 



ds 2 = A 



-dr 2 + {dz - A(s, T)dT) 2 + B(s, T)T 2 dtf + C(s, r)dx 2 
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VI. D7-BRANE EMBEDDINGS: THE DYNAMIC CASE 



To describe the dynamical passage of a cooling quark-gluon plasma through the chiral 
phase transition point, we will embed a D7-brane into the geometry described by the metric 
(I35p . Simulating the evolution of D7-branes in this background is not an easy task. As 
will be seen below, the equations of motion are nonlinear with nontrivial dependence on 
coordinates. Simulations of realistic scenarios will require a serious research in numerical 
methods and will be attempted elsewhere. Here, we consider a simplified scenario and 
test to see if the model at hand gives sensible results, both theoretically and numerically. 
Specifically, we show that a (numerical) solution exists in which the D7-brane pulls out of 
the horizon, indicating the chiral phase transition. We will discuss our numerical methods 
and simplifications in detail, hoping that it will benefit future attempts. 

Like the AdS black hole case, the D7-brane is assumed to fill the five dimensional geometry 
and wrap an 5*3 inside the S5. The brane ends when S3 shrinks to zero size. The three angular 
coordinates of S3 and the five coordinates of the metric (1351) are chosen to be the parameters 
that describe the D7-brane embedding. Throughout this section, we will call the time-like 
direction r the time. Using the previously given form of the five-sphere metric (j2J), we 
assume = and 9 = 6{r,z). This does not take into account anisotropies formed in the 
dual plasma, in a more realistic simulation one would consider dependence on x± and r\. 
The asymptotic behavior of the 9(t,z) field as z — > will be interpreted in the AdS/CFT 
sense. The quark mass M and chiral condensate C of the dual plasma will again be given 
by asymptotic behavior of the scalar field, but now they will be time dependent. 

We set eo = 4/3. The induced metric on the D7-brane is 



ds 2 =\ f 

z l 



(1 - r-lz 4 - z 2 9 2 )dr 2 + 2(z 2 99' - r-h 2 )drdz + (1 + z 2 9 ,2 )dz 2 + r 2 d V 2 + dx\ 
+ cos 2 0d£f. (39) 

The D7-brane action, up to some normalization factor is, 



S D7 = J dzdr J cos 3 9^1 + z 2 9' 2 (l - t'Iz^ - z 2 9 2 - 2r-h 4 9'9. (40) 
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The equation of motion, after some simplifications, becomes 

3r 3 z 2 9 + 3r 3 z A 99' 2 + (3r 5 /V - 3rV) 9" + 3r 3 z A 9 2 9" + (3r 2 z 2 - 3r 7 / 3 z 3 ) 9 
+ (3t 5 / 3 z 5 + r A ' 3 z A + 9t 3 z) 9' + (6r 7 /V - 3rV) 9 3 

+ (Gt^z 11 - 3z 10 - 18r 5 /V + t^z 6 + 12t 3 z 3 ) 9' 3 + (18t 5 /V - 9t a ' 3 z q - 12rV) 9 2 9' 
+ (I8r2 9 - 9t 2 / 3 z 8 - 3Qt 7 / 3 z 5 + 3rV) 66' 2 + Qt 7 / 3 z 4 6' - 6r 3 z 4 96'9' 

- 9t 3 tan 9 + 9t 3 z 2 6 2 tan 6 + I8r 7/3 z 4 66' tan + (9r 5/3 z 6 - 9r 3 z 2 ) 9' 2 tan = 0. (41) 

This is a second order nonlinear partial differential equation. Polynomials of derivatives to 
third order appear, the tangent of 9 contains polynomials of 9 to all orders. The coefficients 
are both space and time dependent. We assume that the Cauchy problem is well-defined 
for this equation, in the sense that for some set of initial conditions defined on a space-like 
hypersurface, the equation will have a unique solution for the causal future. One should not 
expect this statement to hold for every initial condition, we propose consistency requirements 
on the initial condition below. 

In the considered scenario, an embedding of the D7-brane is specified on an initial time 
slice. This embedding is chosen to be an infalling one, corresponding to the high temperature, 
chiral symmetric phase. As the static D7-brane embedding in Lorentzian signature AdS 
black hole background, considered D7-brane embeddings will end at some point inside the 
horizon. Furthermore, boundary conditions are set at the boundary of the space at z — 0, 
the necessity of this will be apparent below. Then the embedding is iterated in time, during 
which the plasma is cooling. Geometrically, this corresponds to the horizon moving away 
from the boundary. We expect to see the D7-brane embedding to be "pulling out" of the 
horizon, through the mechanism discussed above, and at some point becoming non-infalling. 
This will be interpreted as the dynamical passage through the chiral phase transition. 

We start by writing equation (14 ip in first order form. After many trials, we found the 
following form to be numerically more stable. First we define the variables, 

u = 9, v = 9\ (42) 
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and coefficient functions 

a = 3r 3 z 2 + 3r 3 z 4 v 2 , 

b = 6r 7/ V - 6t 3 z 4 uv, 

c = 3r 5 /V - 3rV + 3rW, 

d = (3rV - 3r 7/ V) u + (3r 5/3 z 5 + r 4/3 z 4 + 9r 3 z) v 

+ (I8r^ 9 - 9r 2 / 3 / - 36r 7 / 3 / + 3rV) w 2 + (l8r 5 / 3 z 7 - 9r 4 / 3 z 6 - 12rV) U 2 u 

- 9r 3 tan + 9r 3 ,2 V tan 9 + 18t 7/3 2 4 mw tan 9 + (9r 5/3 z 6 - 9t 3 z 2 ) v 2 tan 9. (43) 

Then, equation of motion ( 14"T|) can be written as 








—u/v 










0' 






u 


+ 


d/ {ya) 


6/a 


c/a 




w' 


= 0. 


(44) 


V 







-1 







v' 







Written in this form the equation is quasi-linear with no source term. The first component 
equation merely states a choice of integrating 9, it actually is a constraint equation. Through 
trial and error we found the above form to be numerically more stable. The second and the 
third component equations contain the real information of the equation of motion, they 
relate the second derivatives of 9. To classify this system as an hyperbolic, an elliptic or a 
parabolic system of equations, one has to solve for the eigenvalues and eigenvectors of the 
3-by-3 matrix appearing on the left of equation HU (lOj. This matrix depends on the solution 
itself, as well as independent variables, therefore it is not possible to do a classification before 
obtaining a solution, which is crucial in posing the correct problem and choosing the numer- 
ical method. On the other hand, the physics of the problem asks for an hyperbolic equation, 
in which the information of the initial time embedding determines the embedding in later 
times (boundary conditions will also be needed as will be discussed below). Therefore, we 
assume that the equation is hyperbolic. This choice puts strict consistency conditions to 
the numerical solution. At every point in space and time, the 3-by-3 matrix appearing on 
the left of equation (|44p should have real eigenvalues (characteristic velocities) with linearly 
independent eigenvectors. The physics of the problem tightens this restriction. Within the 
horizon, we expect the information flow to be only in positive z direction pointing towards 
the singularity. Going back to the equation f T44|) . we see that one of the eigenvalues of the 
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3-by-3 matrix appearing on the left will always be given by —u/v, corresponding to the 
constraint equation. The remaining two eigenvalues (b/2a ± \Jb 2 — Aac/2a) come form the 
second and the third equations, which relate the second derivatives of 9. These eigenvalues 
then must be real and inside the horizon they must be positive. In terms of the variables 
defined above, the reality condition for eigenvalues is (for real 8, u and v ) 

b 2 - Aac> -zV - 2rh 4 uv + rt(l - u 2 z 2 + v 2 z 2 ) > 0. (45) 

For the positivity constraint, the additional condition is 

ac > =► z A + rl (-1 + z 2 u 2 ) > 0. (46) 

In simplifying these equations, we made use of the positivity of r and z. Recalling that the 
horizon is located at Zh(r) = r 1 / 3 , the equation f j46|) is automatically satisfied inside the 
horizon for any embedding. This confirms our expectations. These restrictions should be 
taken into account when choosing an initial condition. 

We use a simple first order upwind scheme to integrate the set of equations fj4*4"j) . Upwind 
scheme is introduced in the appendix |A] 

Having specified the scheme, we now discuss the initial and boundary conditions required 
for solving the system of equations (1441) . Aside from specifying a D7-brane embedding at 
some initial time, the scheme (]A3j) requires boundary conditions at the boundaries of the 
domain to propagate information. This is not just a numerical necessity, the boundary 
conditions carry information from the causal past of the solution that is not included in 
the initial condition. Recall that we consider the scenario in which the initial condition 
is an infalling D7-brane embedding. Initially the D7-brane extends from the boundary of 
the space given by the hypersurface z = to its endpoint z cn< i which is inside the horizon. 
One needs to set boundary conditions at the spatial boundary, z — 0, to accommodate for 
incoming information. The functions 9, u and v will be set at z = and this will be fed 
to the backward differenced term in the upwind scheme flA3j) at the first grid point to be 
propagated in the positive direction. For our numerical simulation, the boundary condition 
at z = is chosen to be 

9(t, 0) = 0, u(t, 0) = 0, v(t, 0) = Mt} /3 , (47) 

where M is a constant quark mass up to some normalization. Note that the temperature 
goes like T ~ r -1 / 3 , hence the r 1 / 3 factor is included in the statement above. 
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At the other hand of the brane, we do not expect a need for a boundary condition. In the 
expected evolution, the endpoint of the brane moves towards the boundary of space, "pulling 
out" of the horizon. The location of the endpoint changes over time, therefore the domain 
of the solution that we are looking for is dynamically being updated at every time step. We 
expect the endpoint of the brane to be determined only by the brane configuration in the 
previous time step. For infalling embeddings, due to our expectancy of characteristic speeds 
being positive inside the horizon there is no numerical necessity for a boundary condition 
at the endpoint. For a non-infalling embedding, we chose to ignore negative characteristic 
speeds at the endpoint and only propagate positive characteristic speeds. We suggest an 
alternative approach below, see footnote ^. Even though one does not expect the necessity 
of a boundary condition at the endpoint, there is a natural consistency condition that comes 
from the equation of motion (HTl) . The brane ends when S3 shrinks to zero size. An inspection 
of the induced metric on the brane (|39|) shows that this happens when 9 = tt/2. The point 
at which the brane ends evolves in time, z en< i = z CQ d(r). However, 6>(r, z en d{r)) = tt/2. Then 



dr 



— ^end^end + #cnd — 0, (48) 

2 = 2 cn d(T") 



where 
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9' " 



i ^cnd r\ 

Z=z end (r) ' ' 



(49) 

2=2end(l-) 



For an infalling embedding, we expect 9 and its first and second derivatives to be nonsingular 
at the endpoint. The most singular terms of the equation of motion fj4TT) are those that are 
proportional to tan 9. Requiring the coefficients of the most singular terms to cancel leads 
to 

1 _ ^2 A2 , -4/3 J3 n' 2 j_ o^-2/3 4 n n> Ji /p- n \ 

1 _ z endPend ' z end° end ~>~ zr z endPendP end z endP end- \ ov ) 

Combining these equations, one can solve for z en d(r), 



z end — T 2 ^ Z end ^ <l/l + ^2 2 ' 
Y ^ end Z end 

We use equations ( )50l) and ( loTT) to check the consistency of our numerical solution. For a 
non-infalling embedding 9 is expected to be singular, so that the brane closes smoothly. To 
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get a constraint on the endpoint one needs to know the exact asymptotic behavior of the 
embedding 8 . 

We choose an infalling embedding as an initial condition. This choice is subject to 
consistency conditions. First of all, it should be compatible with the boundary conditions. 
Moreover, it should obey the restrictions on characteristic velocities discussed above. The 
eigenvalues of the 3-by-3 matrix appearing on the left of equation (jUJ) should have real 
eigenvalues with linearly independent eigenvectors at every point on the grid, which requires 
the condition ( H5|) . The condition of two of the three eigenvalues corresponding to second 
and third equations of the system (JUJ) being positive inside the horizon is automatically 
satisfied, as discussed below equation f )46p . The asymptotics of the embedding near the 
endpoint should satisfy f joTTj) . It is not trivial to find an initial condition that satisfies all these 
properties. Our starting point was noticing that when u = 0, the constraint equation ( 150]) 
on the brane endpoint is satisfied by an infalling D7-brane embedding in an AdS black hole 

1/3 

background with static horizon at Zh = t { . We chose such an embedding, which sets 9 and 
v at the initial time slice, and u = as initial conditions. For consistency, one should look for 
a D7-brane embedding with quark mass that is equal to the boundary condition for the whole 
evolution (which we chose to be static). This is numerically done by a shooting method. 
Moreover, we numerically verified that the reality of characteristic velocities condition (jl5j) is 
also satisfied for this initial condition choice. Figures [3], H] and [5] show the results for a sample 
simulation. Some additional special techniques were used to stabilize the system, which are 
described below. We were able to numerically integrate the D7-brane evolution into a non- 
infalling configuration. After some point numerical instabilities were not controllable, whose 
initial stages can be seen in the figure. 9 

In the figures, the evolution starts from an infalling configuration at Ti = 100 and pro- 

8 Being inspired by the AdS black hole case, we propose the following form of the solution, for z near z on( j 
for an non-infalling embedding 

9(t, z) ~ | - (z cnd (r) - z) 1/2 <*»(t)Cw(t) - *)" (52) 

Although we have not used this proposal in our numerical simulations, in further studies this form could 
be used to propagate the endpoint by solving the equation of motion order by order for z near z on d- 

9 Another set of initial conditions could be given by setting 9, 9 and 9 to zero at some initial time r,, 
which we call the "static initial condition" (although the embedding appears static only momentarily and 
evolves in time), and solve for the embedding at this constant time slice. The equation of motion (|4ip 
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ceeds to a non-infalling configuration. The horizon crosses the endpoint at r = 107.8003. 
After some point numerical instabilities lead to uncontrollable oscillations and the numerical 
integration fails. The initial stages of the formation of these oscillations can be seen in the 
top figure. The quark mass is set to Mr/ 73 = 0.1942. Grid spacing in the z coordinate 
is Az = 0.05. Numerical integration is done until the last point of the D7-brane that is 
on a grid point. Adaptive time steps are chosen to be 0.01 times the maximum time step 
allowed by the CFL condition (IA5I) . In the top and bottom figures we only present 1 time 
step in every 1000 steps. The embedding lines become denser where the numerical integra- 
tion proceeds by smaller step sizes. This was not sufficient to prevent the instability that 
discretization induced on the system. Unstable oscillations grew quickly in time and failed 
;he numerical integration. To prevent this we introduced a sixth order sponge filter, see 
4j for a discussion of sponge filters. Near z = 0, we used a polynomial fit to smooth out 
the unstable oscillations. To check the accuracy of the numerical solution, we checked the 



gives 



3tan0 rr 4/ V + 3 M T 7 5/3 z 2 



-4/3. 



9' -3tan( 



nl2 



+ 2z(2-Tr 4 / 3 z i )6' 3 + li - 9 ri . (53) 

This is to be contrasted with AdS black hole embedding ([5]). One sees that third and the sixth terms 
on the right hand side are new. Let's note that in the z — > z en( j limit, the static initial condition does 
satisfy the consistency condition (|50l) . This actually is a redundant statement since the static embedding 
is solved from the equation of motion anyway, as equation (|50p also is. To find an infalling solution 
for this equation, we follow the same line of argument as in the AdS black hole case. The horizon is 
at Zh(r) = r 1 / 3 . We choose a value for 9 at the horizon that is less than tt/2, 9a < tt/2. Then, from 
the equation above, one can solve for the condition on 9' requiring the coefhecients of the leading order 
singularity to cancel, 

2 ~ 1/3fl ' 3 - -r-X - Ar^X + 3r- 2/3 tan0 o = 0. (54) 



3z 8 - T^V 



-T: 



3 1 u 3 

For an infalling embedding, we expect to have a positive slope at the horizon. The reality and positivity of 
roots depend on particular r,. As in the AdS black hole case, the equation (1531) is numerically integrated 
both towards the singularity and the boundary starting from the horizon with boundary conditions 9a 
and 9q. We found that not all positive real roots of equation ([54| led to boundary conditions that were 
numerically integrable. We tried this setting for a various 9q and Tj and verified numerically the existence 
of initial condition that satisfies the reality condition of characteristic velocities (j45"T) . The time evolution 
of the corresponding initial conditions also led to topology changing D7-brane evolutions, satisfying the 
condition on the endpoint evolution (|50p within numerical accuracy. 



22 



condition on the endpoint (1501) at every time step at the last grid point for the infalling 
configuration. The absolute value of the difference between the left side and the right side 
was always much below the spatial grid size. The error increased when the grid point got 
away from the actual endpoint. For the infalling configuration, the absolute value of the 
error when averaged over all time steps was 0.0002, which we interpreted as a confirmation 
of the numerical result. 

e(x,z) 




100 o 

z 

FIG. 3: The evolution of a D7-brane with the inital condition discussed in the text. 

VII. CONCLUSIONS 

We have made some initial steps toward describing the non-equilibrium dynamics of a 
first order phase transition akin the chiral phase transition in QCD, in a strongly coupled 
supersymmetric large N c gauge theory, by making use of AdS/CFT duality. The dual 
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FIG. 4: The embedding at different time slices. 

description of the phase transition involves a probe D-brane falling into or 'pulling out' of 
an AdS-black hole. 

One aspect of this work involved an attempt to obtain a Lorentzian signature description 
of the latent heat. The probe D-brane geometry in the infalling case implies a lower bound 
on the energy flux at the boundary of Anti-de-Sitter space (and at the horizon) in crossing 
from the high temperature phase to the low temperature phase. However this bound was not 
stringent enough, being less than the actual latent heat computed in Euclidean signature. 
It would be interesting to resolve this discrepancy, which presumably involves an accounting 
for the back-reaction. 

We have also described our initial attempts at simulating the phase transition in a back- 
ground dual to a boost invariant expanding cooling plasma, akin to that produced in heavy 
ion collisions. Our proposed background is of the Painleve-Gullstrand, or 'river-model' form 
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FIG. 5: The same evolution with 8 direction pointing outside the paper plane. In this figure, the 
red (dashed) line shows the evolution of the black hole horizon. 

which seems to be particularly convenient, giving a description valid across the horizon and 
which is manifestly AdS asymptotically. The numerical simulation of the highly non-linear 
equations of motion in a space-time with a complicated causal structure is non-trivial. We 
have restricted ourselves to the simplest case of supercooling, without introducing spatial 
inhomogeneities. In our simulations, we observed a D7-brane 'pulling out' of the black hole, 
however we encountered instabilities shortly after the transition. Further refinements are 
clearly needed to mollify the numerical instability. 
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Appendix A: Upwind Scheme 

Upwind scheme is a method for solving hyperbolic partial differential equations. In this 
scheme, the partial differential equation is discretized by using differencing biased in the 
direction determined by the sign of the characteristic speeds. The direction of propagation of 
information is taken into account through this mechanism. We refer the reader to literature 



for theoretical results on the upwind scheme. We found the discussions in [4], [1 
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37| to be very useful. Here we briefly describe the method, using the notation of [361 ] . 
Given a system of first order partial differential equations in two independent variables, 

^ + B(t,x,y)^ = f(t,x,y), (Al) 

one first introduces a grid over the domain [0, oo) x [0, oo) with grid points (t n ,Xj) defined 
by 

n— 1 

Xj :=jAx (j = 0,1,2,...), T = ^At fc (n = 0,1,2,...). (A2) 

fc=0 

with Ax the spatial grid size and At k the time step. The time step is adaptive, it will be 
defined at each step by the Courant-Friedrichs-Lewy (CFL) condition as discussed below. 
The matrix B is assumed to have real eigenvalues and a complete set of eigenvectors at every 
point over the domain for the particular y in question; this is the condition of hyperbolicity. 
The upwind scheme is given by 

^(yf 1 - y?) + -k B 7 + W ~ yU) + '~(y^ - y?) = f ^ y?)- ^ 

The scheme is solved iteratively for y" +1 - The matrices B™' + and B™'~ are defined locally 
by 

B™' + := S]A]' + (S])-\ I*;' : S'/Aj' (S?) '. (A4) 



26 



S" is the similarity transformation matrix that consists of the local eigenvectors of B" as its 
columns, i.e. if b™ 1; b™ 2 , and b™ 3 are a complete set of (right) eigenvectors of B" for the grid 
point {t n , Xj) and the value of the function y at the grid point, y™, then S™ = (b™i, b™ 2 , b™ 3 ). 
A"' + and A™'~ are diagonal matrices that contain the positive and negative eigenvalues of 
the matrix B™ as diagonal elements respectively, i.e. if A™ 1; A™ 2 and A™ 3 are the eigenvectors 
corresponding to b" )1( b™ 2 , and b™ 3 , then (A™' + )jj = max(A™j, 0) and (A™'~)jj = mh^A^O), 
i = 1,2,3. B™' + and B™'~ give a decomposition of B™, B™ = B™' + + B™'~. We note again 
that this decomposition has to be done at every grid point and will be different for different 
functions y. 

The eigenvalues A? l3 A™ 2 and A™ 3 are the local characteristic velocities of the system, 
they describe the direction and speed of information flow. The upwind scheme (1A3I) uses 
a backward differencing to propagate information in the positive direction and a forward 
differencing to propagate information in the negative direction. The time stepping should 
be adaptively chosen by the CFL condition, 

At n 

max(A^)-— < 1, (A5) 



Ax 



This condition prevents acausal information flow on the 



for stability (4, 
grid. 

The upwind scheme is dissipative. There are many other schemes with better accuracy, 
but their implementation to the equation at hand is very difficult due to the nonlinearity 
and parameter dependent coefficients. 
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